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Abstract: We apply a dressed perturbation theory to better organize and economize the 
computation of high orders of the 2-body effective action of an inspiralling Post-Newtonian 
gravitating binary. We use the effective field theory approach with the non-relativistic 
field decomposition (NRG fields). For that purpose we develop quite generally the dressing 
theory of a non-linear classical field theory coupled to point-like sources. We introduce 
dressed charges and propagators, but unlike the quantum theory there are no dressed 
bulk vertices. The dressed quantities are found to obey recursive integral equations which 
succinctly encode parts of the diagrammatic expansion, and are the classical version of the 
Schwinger-Dyson equations. Actually, the classical equations are somewhat stronger since 
they involve only finitely many quantities, unlike the quantum theory. Classical diagrams 
are shown to factorize exactly when they contain non-linear world-line vertices, and we 
classify all the possible topologies of irreducible diagrams for low loop numbers. We apply 
the dressing program to our Post-Newtonian case of interest. The dressed charges consist 
of the dressed energy-momentum tensor after a non-relativistic decomposition, and we 
compute all dressed charges (in the harmonic gauge) appearing up to 2PN in the 2-body 
effective action (and more). We determine the irreducible skeleton diagrams up to 3PN 
and we employ the dressed charges to compute several terms beyond 2PN. 
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1. Introduction and Summary 

Gravitational wave observatories (see for example the reviews ]l]] and references therein) 
demand knowledge of the waveform emitted by an inspiralling binary system of compact 
objects. Fully Generally Relativistic numerical simulations can now simulate such wave- 
forms - see the review [Q] and references therein. Yet, as always, an analytic treatment 
is complementary and improves insight, especially into the functional dependence of the 
results on the parameters. A perturbative analytic treatment is possible in two limits. The 
first is the Post-Newtonian (PN) approximation, and it holds whenever the velocities are 
small compared with the speed of light, or equivalently through the virial theorem when- 
ever the separation between the compact objects is much larger than their Schwarzschild 
radii. The second limit is that of an extreme mass ratio. In this paper we shall concentrate 
on the PN approximation which is always valid at the initial stages of any inspiral. 

The computation of the effective 2-body action in PN was a subject of considerable 
research over the last decades and the current state of the art is its determination up to 
order 3.5PN, as summarized in the review || (see also the recent [Q]). Another approach, 
the effective field theory (EFT) approach to this problem was suggested by Goldberger and 
Rothstein (2004) g, where more traditional GR methods are replaced by field theoretic 
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tools including Feynman diagrams, 1 loops and regularization. In particular ||] reproduced 
the 1PN effective action (known as Einstein-Infeld-Hoffmann) within the EFT approach. In 
0, [8| the metric components in the Post-Newtonian, non-relativistic limit were conveniently 
decomposed into a scalar Newtonian potential, a gravito-magnetic 3-vector potential and 
a symmetric 3-tensor. These fields were termed collectively "fields of Non-Relativistic 
Gravitation (NRG fields)", and the derivation of 1PN was shown to further simplify. In 
the 2PN expression was reproduced within the effective field theory approach together 
with NRG fields. 

Other recent developments related to either PN or the EFT approach include: dissi- 
pative effects and effective horizon degrees of freedom |I~0|, 11]; thermodynamics of caged 



black holes 0, [b| through the EFT approach ||, [7|, [l§|; EFT |§ (see also 0, |lj§) 
Hamiltonians [19, 2C] for rotating point-particles; tidal effects for compact objects [[IT], 22]; 



approximate solutions for higher dimensional black object including rings [23] and black- 



folds [24|| ; a mechanized EFT computation for 2PN [25|; "De-Turek" gauge for numerical 



relativity [26]; and finally radiation reaction and waves within the EFT approach |27|, 28]. 

The computation of order 2PN || demonstrates a proliferation in the number of dia- 
grams: from 4 at 1PN to 21 at 2PN, and furthermore the number is expected to continue 
and grow at the next order. Clearly, it would be useful to have an improved perturbative 
expansion. In this paper we shall present a new method to better organize the calculation 
and economize it. The basic idea is to recognize recurring sub-diagrams which physically 
describe "dressed" charges and propagators. The simplest example is furnished by the 
Newtonian potential Gm\m,2l\r\ — T2\ C S e ff which belongs to order OPN. This term is 
associated with the diagram in figure |l|. At higher orders there is a class of diagrams of the 
form shown in figure |2] where the point masses are replaced by dressed energy distributions 
defined on the top row of figure || 

mi5(Xi - Ti) -> p dr (xi - n) (1.1) 

and the bare propagator is replaced by the dressed (relativistic) one defined on the bottom 
of figure |3] 

p S(At) - G(k, At) . (1.2) 

The same sub-diagrams which define the dressed energy distribution and appear in the 
calculation of the dressed Newtonian interaction (figure Q) appear also in other diagrams. 
It makes sense to record the values of these sub-diagrams and later re-use them. This is 
the basic idea of the dressed perturbation theory. 

In addition, the dressing procedure is shown to economize the calculation in a more 
significant way as follows. We find that the dressed couplings satisfy a certain recursive 
integral equation schematically described in figure A perturbative expansion of the 
solution to this equation equals an infinite sum of diagrams, and in this sense it encodes 
many diagrams. This is nothing but the classical version of the Schwinger-Dyson equations 



1 See also M for certain early days efforts to calculate classical potentials to higher post-Newtonian orders 
using field theoretical methods and Feynman graph techniques. 
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Figure 1: The diagram which represents the Newtonian potential interaction mediated through 
the Newtonian scalar field <f>. The notations will be fully defined later in section pu 



Figure 2: A class of diagrams that we interpret as the Newtonian potential interaction between 
dressed energy distributions of each body through a dressed propagator. The dark blobs represent 
any sub-diagram with an arbitrary number of vertices on the world-line and a single external leg 
for the Newtonian potential. There are no bulk loops, as always in classical physics. The light blob 
represents any sub-diagram with two external legs of the Newtonian potential, which amounts to a 
propagator with an arbitrary number of retardation insertions. 



Figure 3: The diagrammatic definition of the dressed energy distribution (top) and the dressed 
propagator (bottom). The dressed energy distribution is defined through the one point function 
for <fi in the presence of a single source (after stripping external propagators), while the dressed 
propagator is defined through the full two point function for <f>. 

which were first written in the context of quantum electro-dynamics [p9f| . Yet the current 
classical version is of a higher practical value compared to its quantum counterpart since 
it does not involve the infinitely many dressed bulk vertices. 

Some discussions of classical versions of the Schwinger-Dyson equations appeared al- 
ready, yet they all appear to consider a significantly different context. studied unequal 
time correlation functions of a non-equilibrium classical field theory, while aims at 
giving a construction of the local algebras of observables in quantum gauge theories. 
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Figure 4: A schematic diagrammatic representation of the Schwinger-Dyson recursive integral 
equation satisfied by the dressed couplings. More details are given in the body of the paper. 



The paper is divided into two parts. In section || we describe and discuss the dressed 
perturbation theory for a general classical (effective) field theory, while in section ||] we 
apply it to the Post-Newtonian theory. We start in section [2] by considering a simple 
scalar classical field theory coupled to point-like particles as the context for introducing the 
required dressing concepts. In subsection [2 . 1| we characterize the fully factorizable diagrams 



before turning in subsection 2.2 to our main definition, that of the dressed quantities. In 
subsection ^3] we define the dressed perturbation theory and assert that it is equivalent to 
the bare theory. 



In subsection 2.4 we explain the recursive integral equation a la Schwinger-Dyson and 



comment on its relation with related concepts. We close the general theory section in 



subsection |2.5| with a classification of irreducible skeletons at low loop numbers - diagrams 
which are neither factorizable nor do they include dressed sub-diagrams. 

In section we apply these concepts of dressing theory to better organize and econ- 



omize the two-body Post-Newtonian (PN) effective action. In subsection 3.1 we set up 
the problem by displaying the Post-Newtonian effective action and the associated Feyn- 
man rules. In subsection |3.2| we explicitly compute the three charges in PN. The charge 
which couples to the gravitational field is the energy-momentum tensor, and since the 
gravitational field naturally decomposes in the PN limit into three NRG-fields the source 
decomposes too into three corresponding parts: the energy density, the momentum density 
and the stress. These charges are computed both in fc-space and in position space up to the 
order required to reproduce the 2PN results, and the stress is determined to one additional 
order. 

In we analyze qualitatively the diagrams relevant to the computation of the two- 
body effective action. Starting with OPN and making our way to order 2PN we explain how 
some of the diagrams can be expressed through dressed charges, and which diagrams can 
be generated by the integral equation. Up to order 2PN only a single (non-tree) diagram is 
found to be irreducible. We confirmed that the known 2PN effective action is reproduced 
after incorporating our expressions for the dressed charges from the previous subsection. 
While thus far almost all the computations are related to the known [|| 2PN effective action 
which was already reproduced within the EFT approach with NRG fields in ||, here we 
list all the skeletons required for the computation of 3PN, thereby indicating the road-map 



for organizing the 3PN computation. In subsection 3.4 we convincingly demonstrate the 
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utility of our method by computing certain novel 3PN and 4PN diagrams. 2 Finally, in 
appendix |A| we collect some useful integrals. 



2. Dressed perturbation theory in CLEFT 

In this section we define a dressed perturbation theory in the general CLEFT (Classical 
Effective Field Theory) context, namely for any classical non-linear field theory coupled to 
point-like sources. 

In order to illustrate the main ideas in a simple setting we consider a scalar field model. 
The generalization to an arbitrary CLEFT is straightforward and will be spelled out at the 
end of subsection |2]^ on p. ||. 

Consider the following bulk action 



SbuikM = d x 



- 2 W +^-p 2 -- s 



(2.1) 



for a scalar field (ft with propagation speed c and coefficient of cubic interaction a, where 
(p couples to any point particle of mass m and charge q through 

S p = -( m - q ) J dr-q J e^ (r)) dr , (2.2) 

dr is the proper time element and the particles will be assumed non-relativistic dx/dr <C c 
(or even static in part of the discussion). The total action for a many body system is 

S = Sbuik + XT ^P' a ' 



where S Pta is the world-line action (2.2) of the a'th particle characterized by m a , q a - 

Let us briefly discuss the considerations for choosing this form of the action. The 
retardation term proportional to 1/c 2 is considered as a small perturbation in the non- 
relativistic limit and represents a general small perturbation of the quadratic term. The 
term proportional to a was chosen to represent any non-linear interaction. The inter- 
action term — q J exp (0(x(r))) includes the charge coupling — q J 4>{x{t)) together with 
some representative non-linear terms (this exponential form of the interaction appears in 
Post-Newtonian theory, for example). 

The bulk theory (2.1) has a vacuum at 4> = and we consider the perturbation theory 



around it. As usual, we note that while this vacuum is unstable, it could be stabilized by 
adding a mass term to /S^uiki an d it could even be made a global minimum through the 
addition of a quartic term in the potential. 

In this paper we study the two-body problem rather than the seemingly more general 
many body problem since it is the simplest and currently the most interesting case. The 
generalization to the many body problem seems straight forward. 

The Feynman rules involving (f> are shown in figure [|. They are standard except for 
the CLEFT conventions which makes them real - vertices are read from the action and 
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Figure 5: The Feynman rules involving <j> for the static scalar field theory whose action is given 
by ( |2.3| ). The rules describe: the propagator and the quadratic perturbation vertex - retardation 
(top), the cubic bulk vertex (middle), and the world- line vertices (bottom) where the ellipsis stand 
for additional non-linear world-line vertices. Hereafter k denotes a spatial wave-number. 

propagators are given schematically by —I/S2 where S2 is the quadratic part of the action. 
See [^] for the original definition and full detail. 
The two body effective action is defined by 

S e ff[x 1 ,x 2 ] := S[x 1 ,X2,4>(xi,X2)\ (2.4) 

where the RHS represents the action ( |2.3[) evaluated on the solution (j) = 4>{x\,X2) given 
the particle trajectories x\{j\), X2{t2)- The effective action is known to be equal to the 
sum of all connected Feynman diagrams made out of (f> propagators with arbitrary bulk 
vertices and world-line vertices but without loops of propagating fields (such classically 
forbidden loops are allowed quantum mechanically). 

2.1 Factorizable diagrams 

There are several possible paths to classical dressed perturbation theory. We shall build it 
from first principles and later discuss its relations with both the quantum version and the 
standard classical theory. 

The main idea is to economize the perturbation theory by identifying certain recurring 
sub-diagrams. Before we proceed to the more essential, dressed sub-diagrams we discuss a 
stronger and simpler form of reduction, namely factorizable diagrams. 

A Feynman diagram is called factorizable whenever the expression which it represents 
factorizes (into a product of factors), each one corresponding to a sub-diagram. An example 

2 The computed first Post-Minkowskian approximation |^2| includes information about a certain class of 
diagrams to arbitrarily high PN orders, which are different however from the ones we compute here. 
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Figure 6: An example of a factorizable diagram (the simplest). The two factors are circled (in 
red) and they intersect in a NL-WL vertex. A diagram of this type appears at 1PN - see figure |2^. 



is shown in figure |(| In CLEFT we have the interesting property that a diagram of the 
2-body effective action is factorizable if and only if it contains a Non-linear world-line 
(NL-WL) vertex, where a NL-WL vertex is a vertex with more than a single bulk field. 

Indeed, if a diagram contains a NL-WL vertex then since quantum loops are not 
allowed, each leg of the vertex generates a separate sub-diagram. Conversely, a connected 3 
factorizable diagram necessarily factorizes at a vertex. A bulk vertex would not serve since 
wave-number conservation couples between all its legs. 4 

2.2 Dressed charge and propagator 

Following our discussion of factorization we proceed to consider only diagrams without any 
NL-WL vertices. Any non-factorizable diagram contains various sub-diagrams. We wish 
to further economize the perturbation theory by identifying in a unique and natural way 
a class of sub-diagrams which repeatedly show up at high orders of the two body effective 
action. We shall call them "dressed sub-diagrams". We start with constructive definitions 
and an explanation of the name's origin, to be followed by a more abstract characterization 
through an equivalence of perturbation theories which serves to explain the rational behind 
the definitions. 

Definitions 

• The dressed charge of the particle, Pdr(k, t) is defined diagrammatically through figure 
(top). The dressed charge is an infinite sum of sub-diagrams, where each summand 
will be called a dressed charge sub-diagram. 

• The dressed propagator, Gd r (k,At) is defined diagrammatically through figure 
(bottom). Each summand in the definition will be called a dressed propagator sub- 
diagram. 

3 AU the diagrams of the 2-body action are connected by definition. 

4 In CLEFT we cannot have a factorizable diagram such as the "figure 8" diagram in (j> 4 theory, since we 
cannot have k conservation for any strict subset of propagators leaving the bulk vertex, as each propagator 
connects to a distinct world-line vertex (or vertices) where k is arbitrary. 
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Figure 7: The diagrammatic definition of the dressed charge distribution (top) as the one point 
function for <p in the presence of a single source and the dressed propagator (bottom) as the full 
two point function for (j) in vacuum. Retardation vertices are forbidden on the external leg of the 
dressed charge. 



Figure 8: The diagrammatic representation for the value of the field in the presence of the point- 
particle source. The only difference with the definition of the dressed charge in figure [?] (top) is an 
added propagator on the external leg. 

Let us inspect these definitions. In equations Pdr{ r ) is denned through 



where the second equality is a differential equation which together with retarded boundary 
conditions defines <3?(r, t), the full solution for the field eft in presence of the given source 
world-line. This equation is equivalent to the diagrammatic definition since <3? is given 
diagrammatically by figure ^, and hence in /c-space &(k,t) = —pd r {k } t) /k 2 from which 
J2.5D follows in position space. 

The dressed charge describes the apparent particle charge (distribution) at long dis- 
tances. It arises from the non-linear interactions of the scalar field which "dress" the point 
charge, and is useful for studying the dynamics of a system composed of several such par- 
ticles. The term which we use, "dressed" charge (and propagator) is standard terminology 
in quantum field theory (see for example [p3f| ) as well as in classical field theory. For com- 
pleteness we would like to mention some related terms. The first is renormalization, where 
like here one replaces a bare quantity by a scale-dependent renormalized quantity which 





(2.5) 



S 



is defined through the divergences of the same sub-diagrams. 5 Another related term is 
resummation where one discusses partial sums of diagrams. 
The equation form of the dressed propagator is 

r 1 - 1 - 1 1 

dr □ k 2 + d 2 /c 2 k 2 i + d 2 /{c 2 k 2 ) 

-f( 1 -3P* + ?F* + -) (2 ' 6) 

where the last equality is a representation of the series on the left of figure ^ (bottom), 
while the previous expressions can be considered to be a closed form summation of that 
series. 

Physically, in this field theory and also in PN the dressed propagator is nothing but 
the fully relativistic propagator. However, this need not be the case in general. 



Example. To illustrate these ideas let us discuss an example. A detailed application to 
the Post-Newtonian theory will be given in section ^. 

Consider the 6-loop Feynman diagram in figure |9|(a) which contributes in the bare 
theory to the 2-body effective action. It has five dressed sub-diagrams (circled) - two 
dressed charges and three dressed propagators. Note that all the world-line vertices of 
each dressed charge belong to one and the same point particle. Note also that retardation 
vertices are allowed inside the dressed charge but not on its external leg. Upon replacing 
the dressed sub-diagrams by dressed vertices and propagators we obtain the diagram's 
skeleton 6 in (b). The resulting diagram is only two-loop. The other loops were absorbed 
by the dressed charges. Note that the 3-loop dressed charge sub-diagram includes in it 
other dressed sub-diagrams, but these are not maximal. 

Generalization. So far we worked in a simple setting of a cubic scalar field model. Here 
we shall indicate how the definitions given above for the dressed sub-diagrams generalize 
to a general field theory. 

In the scalar field theory we had a single dressed vertex and a single dressed propagator. 
In a general Classical field theory interacting with point-like particles we have a dressed 
charge for each field. For the dressed charge to be non-trivial the bulk theory must be 
interacting. 

The dressed propagator in a general field theory is labeled by any two fields that can 
"mix". In case the fields do not mix then we have a single dressed propagator for each 
field. For the dressed propagator to be non-trivial we need the quadratic Lagrangian to 
decompose into a leading part and a perturbation, so that the leading part will determine 

5 The idea of renormalization in the PN context appeared already in the literature. associated 
renormalization with the regularization of certain divergences which appear in the particle's effective action 
at order 3PN. |eJ studied the same divergences from the EFT approach and in particular they computed the 
renormalization of the energy-momentum tensor of a static particle up to 2 loops. Later |}5| also studied 
the renormalization of the energy-momentum tensor. 

6 This term will be further discussed and defined in the next subsection. 
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(a) (b) 

Figure 9: An example for the correspondence between diagrams of the bare theory, such as (a) 
and the corresponding dressed diagram (b). Note a standard subtlety: retardation insertions are 
allowed inside a dressed charge vertex, but not on its external leg. See text for further discussion. 



the bare propagator, while the small part will determine the two-vertex. The dressed 
propagator will then be proportional to the inverse of the full quadratic Lagrangian. We 
note that one may choose to diagonalize the dressed propagator and accordingly redefine 
the fields such that there will be no mixing through dressed propagators. 

2.3 Equivalence of perturbation theories 

We proceed to define a dressed perturbation theory, which is equivalent to the original one, 
but somewhat more economic. 



Definition. A Feynman diagram which includes a non-trivial dressed sub-diagram of the 
form shown in figure |^ (dressing sub- diagram) will be called dressing-reducible. Otherwise 
it will be called dressing-irreducible. 



Definition. The dressed perturbation theory is defined as follows 



• Figure 10 shows the changes in Feynman rules relative to the original theory (figure 
I)- 

• Only dressing- irreducible diagrams are allowed. 

The original perturbation theory will be distinguished from the dressed one by referring 
to it as bare. 



Property - The bare and dressed perturbation theories are equivalent: each diagram of 
the bare theory is included exactly once in a dressed diagram. 

This property is analogous to a standard one holding for dressed actions in QFT, and 
we shall outline a proof. Given a bare diagram we claim that its maximal dressed sub- 
diagrams are unique. We find the uniqueness property to be quite apparent when one 
thinks about it, but we shall not attempt to provide here a proof, as it seems tedious. 
Given the decomposition we replace the dressed sub-diagrams by propagators and vertices 
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-Q -p dr (t,k) 




Figure 10: The dressed Feynman rules contain changes relative to the bare Feynman rules in 
figure ^. The bare world-line vertex (charge) Feynman rule changes into the dressed one, the bare 
propagator changes into the dressed one, and the retardation rule is omitted. The rest of the rules 
remain unchanged. 



of the dressed theory according to figure 10. The resulting reduced diagram is called the 
skeleton of the original diagram. The uniqueness of decomposition implies now that each 
bare diagram is included once and only once in the dressed perturbation theory, and hence 
the two are equivalent. 



Discussion. Before proceeding to describe another property of the dressed theory, namely, 
the integral equation, let us pause to discuss some aspects of the definitions and property 
above. 

Rational behind definitions. The decomposition into a skeleton with blobs which repre- 
sent dressed sub-diagrams is natural in a practical, computational sense. When computing 
a diagram, the dressed sub-diagrams are almost inevitably evaluated on the way, by their 
nature. Hence it is computationally natural to prepare a list of dressed sub-diagrams and 
their value, in order to avoid their repeated evaluation, which is exactly what the dressed 
theory does. 

Yet, the computational argument alone does not fix our definitions of the dressed sub- 
diagram. The characterizing property is precisely the equivalence property above, namely 
that the dressed perturbation theory is equivalent to the bare one, or in other words, the 
dressed theory is self-sufficient or autonomous. We claim (again without proof) that this 
property can be used to derive our definitions. 

Analogy with effective action in Quantum Field Theory ( QFT). The ideas above are 
analogous to those of the effective action in standard QFT. There one collects all the 1- 
particle- irreducible (1PI) diagrams into the effective action and then allows in the perturba- 
tive expansion only tree diagrams of the effective action. In both cases it is important that 
the decomposition is unique - any diagram can be uniquely divided into 1PI sub-diagrams 
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Figure 11: An over-counting problem with an alternative definition of the (non-static) dressed 
perturbation theory, as discussed in the text. 

which allow a reduction to a tree skeleton. 

Relation of the dressed sub-diagrams with the "dressed source" in standard QFT. The 
two notions are essentially the same. Note however that unlike some cases, for us it is 
important that the dressed charge sub-diagrams are only those where all world-line vertices 
belong to one and the same point particle. 

Our notion of the dressed propagator is also essentially the same as the dressed prop- 
agator and the associated field strength renormalization in standard QFT. The 2-vertex 
in CLEFT is the "self energy" (the 1PI two point function), and (|2.6| ) is essentially the 
standard QFT relation between the self energy and the dressed propagator. The difference 
is that in QFT many diagrams can contribute to the self energy (actually normally they are 
infinitely many corresponding to an arbitrary number of possible loops) while in CLEFT 
the 2-point vertex is read directly from the Lagrangian. 

Why is the dressed propagator necessary? Consider doing away with the definition of 
the dressed propagator by allowing the dressed vertex sub-diagrams to include retardation 
vertices on the external leg. In the classical theory this would actually work in all but one 
important class of diagrams, that of the dressed Newtonian interaction - figure |2], where 
the correspondence between the bare and dressed diagrams would break. For example, 



consider the diagram in figure |11|. There are two distinct ways to cut the diagram into sub- 
diagrams, the cuts being denoted by a,b. Accordingly this diagram is doubly counted in 
the putative dressed theory, thereby disqualifying it. This is all the better since the dressed 
propagator is an appealing physical concept which we would not want to lose anyway. 

2.4 The Schwinger-Dyson recursive integral equation in CLEFT 

In this section we shall describe a second property of the dressed perturbation theory: 
certain recursive relations which generally take the form of integral equations. We start by 
considering the static limit of the scalar theory fl2.3| ) where the idea is simpler to illustrate 
and later we refine it to include the general non-static case. 
Recall the definition of pdr(f) in figure |7[ 

Property. The dressed quantities satisfy recursive relations. In the static limit pd r ( r ) 
satisfies the recursive relation which is shown diagrammatically in figure 12. Its equation 



- 12 - 



Figure 12: The diagrammatic representation of the recursive integral equation satisfied by the 
dressed charge in the static limit. 



form reads 



d 3 k a f d 3 k f d 3 fa p dr (fa) Pdr{k - fa),, 

- K pdr{k)<p(-k) = - q<j)\? =0 - - I j^—^ I j^—^ — u 71 rvr<P{- k ) ■ 



(27r) 3WVVrv ; ^ lr=0 2 J (2vr) 3 J (2vr) 3 k\ (k - faf 

(2-7) 

After factoring out J d 3 k (j)(—k)/ (2tt) 3 we are left with the following integral equation for 

Pdr(k) 

f t \ « f d 3 fa p dr {fa) p dr (k - fa) 

PdA.k)-q+-j — (fc-fci) 2 ' { } 



Given a small a the integral equation can be solved perturbatively in a by expanding 
Pdr(k) = a n p( n \k). The zeroth order is given by p(°\k) = q, the first order is given 
by p^\k) = i J d 3 fa/(2ir) 3 a n^^yi and so on. This iterative solution of the integral 
equation is precisely equivalent to the diagrammatic expansion in figure |7| (because they 
both compute the same quantities, namely p^ n \k)). 

The advantage of the integral equation over the diagrammatic expansion is that it is 
shorter and economizes the computation by avoiding the need to identify all the necessary 
diagrams and compute them. Some readers may benefit from the following analogy: the 
relation between the recursive relation and the full diagrammatic expansion is analogous 
to the relation between defining a function through a differential equation and defining it 
through the corresponding power series which solves the differential equation. 

The non-static case. Once we restore time dependence into the action ( p.3| ) we expect 
the dressed propagator to have a role in the recursive relations as well. First, it has a 
recursive relation of its own given diagrammatically by figure |l3| (bottom). In this case 
the recursive relation is actually algebraic rather than integral, and hence the closed form 
solution Q2.6| ) exists. Secondly, the recursive relation for the dressed vertex (figure |l2| ) is 
refined to include the dressed propagator and is now given by figure |l3| (top) . 

We now proceed to make several comments. 

• The recursive relation (figures [l^|l3| and equation |2.8| ) can be considered as inherited 
from its quantum version (see for example |36[| ), but there are interesting differences. 
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Figure 13: A diagrammatic representation of the full set of (the two) recursive equations satisfied 
by the dressed charge (top) and propagator (bottom) in the general, non-static case. Note a change 
in the recursive relation for the dressed charge relative to the static equation in figure O. 



As it is, figure 12 does not hold in quantum field theory 7 because there can be 



additional 4> propagators connecting the two blobs on the RHS. Since such additional 
propagators would create a closed loop of propagating fields it is forbidden in CLEFT. 
In order for the recursive integral relations to apply in the quantum case one must 
generalize them to include the dressed bulk vertices yielding the celebrated Schwinger- 



Dyson equations 1 29 ] . However, as stated in several textbooks the practical usefulness 
of the quantum version is rather limited since more and more Green's functions with 
more external legs participate in the equations as the order is increased 8 while 
in CLEFT the equations are more practical exactly because there are only a finite 
number of dressed quantities. 

The dressed mass pdr(f) is closely related to the field profile <fr(r) generated by the 
point-like source. Indeed, each one of them is sufficient to determine the other through 
( |2.5D . While pdr( r ) satisfies an integral equation, <5(r) satisfies a "mirror" differential 
equation, namely the equation of motion given by the second equality in ([^5|). The 
perturbative expansion of the integral equation is dual in turn to the perturbative 
expansion of the differential equation into some sort of a power series (which will 
generally include log factors as well). 

Relation with the beta function. The recursive integral equation (|2.8|, and figures 



7 except for the tree diagram limit of course. 

8 "From the point of view of making a practical calculation we have accomplished little; the unknown 
quantities ... have been expressed in terms of yet another unknown..." [Q, "... the system involves an 
infinite hierarchy of equations ... their usefulness is limited" {^J. 
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Figure 14: The top line shows a possible definition for a 2-field dressed world-line vertex. The 
second line demonstrates the over-counting problem which occurs. The bare diagram on the left 
can be cut into sub-diagrams in two different ways denoted (a) and (b). 

[T^]l3|) determines Pdr{k) and so does the beta function. Yet, the two equations are 
different as the beta function is a first order differential equation for dpd r / d\og(k). 
Therefore it must be that the beta function equation is a special or limiting case 
(whose precise definition will not be pursued here) when the leading behavior of pd r 
is logarithmic in k. 



2.5 Irreducible 2-body skeletons 

Consider the non-linear world- line vertices (NL-WL vertices) in the scalar action fl2.3|), 
namely vertices with more than a single bulk field. Our definition of the dressed charge 
concerns a world-line vertex with a single bulk field. It is possible to generalize the definition 
of the dressed vertex to any number of bulk fields, such as the 2-fleld vertex shown on the 



top line of figure 14. In this way one may define a fully dressed one body effective action. 

However, for the purpose of computing the 2-body effective action we use only the 
dressed charges (and propagators) and for the NL-WL vertices we use the bare vertices 
rather than the the dressed ones because they would have created a problem: the decom- 
position of diagrams would not have remained unique (just as in our discussion around 
figure 11). For example, the diagram on the second row of figure |14| can be decomposed 
in two different ways by the two shown cuts, and that corresponds to doubly counting this 
diagram in the putative fully dressed perturbation theory. 

Recall (subsection |2.1| ) that non-linear world-line vertices are interesting for another, 
complementary property: a diagram of the 2-body effective action is factorizable if and only 
if it contains a Non-linear world-line vertex. From this perspective the above-mentioned 
issue with using dressed WL-NL vertices is all the better since there is no need for dressing 
- all diagrams with NL-WL vertices are factorizable and hence reduce to computations of 
lower order and lower loop number. 



2-loop : 



3-loop : 




Figure 15: Classification of topologies for irreducible diagrams of the 2-body effective action, both 
non-factorizable and dressing-irreducible. The top line shows the possibilities at 2-loop, and the 
bottom line shows the possibilities at 3-loop. 



It is interesting to classify the possible irreducible diagrams in the 2-body effective action 
- those which are both non-factorizable and dressing-irreducible. The possible topologies 
are independent of the details of the specific field theory. At 1-loop there are no irreducible 



diagram topologies. At 2-loop (top line of figure 15) there is the "H" diagram on the left 



and its degeneration - the "X" diagram where the two cubic bulk vertices degenerate into 



a quartic vertex. At 3-loop (bottom line of figure 15) there is the topology shown on the 
left, together with its two possible degenerations. We did not list here the trivial 0-loop 
(tree-level) diagram shown in figure |l[ 

3. Dressing the 2-body Post-Newtonian problem 

In this section we apply the dressed perturbation theory to the Post-Newtonian (PN) 



expansion, thereby demonstrating its utility. We start in subsection by setting up the 
problem, establishing the conventions, and displaying the Post-Newtonian effective action 
and the associated Feynman rules. In subsection |3.2j we explicitly compute the three PN 
charges both in A:-space and in position space up to the order required to reproduce the 2PN 
results, apart from the stress that is determined to one additional order. In [T^ we analyze 
qualitatively the diagrams relevant to the computation of the two-body effective action 



including a list of all the skeletons required for the computation of 3PN. In subsection 3^ 
we compute certain novel 3PN and 4PN diagrams. 

3.1 Effective action and Feynman rules 

Consider a binary system composed of slowly moving massive objects. Replacing the 
vicinity of each of these objects by the relativistic point particle coupled to gravity and 
neglecting the effects of the higher-dimensional tidal terms, leads to the following effective 
action 

S = Seh + S pp . (3.1) 
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In PN the gravitational metric field is naturally divided into three fields of Non- 
Relativistic Gravity (NRG-fields): (ft the Newtonian potential, Ai the gravito- magnetic 
vector potential and 7^ a symmetric tensor (with spatial indices) |^. The field redefinition 
9tw -> ((ft, Aimj) is defined by @, |, § 

ds 2 = dx" dx v = e 2<t> (dt - Ai dx 1 ) 2 - e -W~ 3 ) 7ij dx i dx j , (3.2) 

where we take the space-time dimension d to be arbitrary and only at the end we shall 
specialize to 4d. The bulk action is the Einstein-Hilbert action expressed in terms of 
NRG-fields 

s EH [g) = ^ J R[g] - s EH [ 7 , a, 4>] . (3.3) 

The point particle trajectory is denoted by x = x(t) and we shall denote by v = x(t) 
and a = x(t) its 3-velocity and 3-acceleration, respectively. The point particle action is 
given by 

S PP = "E ma / dT = -Y j m a f dt S> y/(l-A- v a ) 2 - e -2(d-2)*/(d-3) ~ w i v { = 

a J a 

= ~ ma j dt y /l ~ V a 



^ I \(d-3)^r^J yr^i 271^2 1 

where the dummy index a runs over all the masses involved in the binary evolution, in the 
second equality equation (3.2) was applied, the ellipsis denote terms non-linear in the bulk 
fields ((ft, A, 7) and we define 

&ij ■= Hj ~ Sij . (3.5) 

The Einstein-Hilbert action is invariant under reparametrizations and should be gauge 
fixed. Leaving aside the question which gauge would be optimal for this problem, we choose 
the harmonic gauge in order to facilitate comparison with the literature. Accordingly, we 
add the following gauge fixing term to the Einstein-Hilbert action 

Sgf = 3^ J d d x^g- , (3.6) 

where = T^g^ . 

The Feynman rules for (ft, Ai and Uij coupled to the world-line can be read from fl3.^ ) 
and we list on figure 16 those couplings which are necessary to our discussion. Solid, 
dashed and wavy lines of the figure are associated with propagators of the instantaneous 
non-relativistic modes (ft, Ai and <Jij, respectively. In momentum space these propagators 
are given by 

= 8irG ^Ef S(t - ^ 1 



= -WttG 5(t - t') 5ij £ 

VWV = Z2irG5(t-t')P lhkl ^ 
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171 J ale {d-3)y/l^ 



Figure 16: Feynman rules obtained from the expansion of (jM) up to linear order in tfi, Ai and 



u i j . The undetermined wave- numbers flow into the vertex. 



with Pj 



ij,kl 



Sik$jl + $u5jk 



d-3°W 



Ski 



The Feynman rules for the bulk vertices are obtained from expansion of the Einstein- 
Hilbert action (|3.3|) and the gauge fixing term (|3.6|). The resulting set of vertices which 



contribute to the calculations below are presented in figure 17 and figure 18, where we 



separate vertices which involve time derivatives in figure |18j, from the static ones in figure 
17. 

In addition, one has to impose wave-number 9 conservation at each bulk vertex by 
assigning it a delta-function factor, (2vr) d_1 (5(^,- kj), where kj is the total wave-number 
flow into a given vertex, and one must integrate over each undermined wave-number k of 
the diagram 

d d - l k 



L-S 



\d-l 



Finally, one has to divide by the symmetry factor of the diagram. 



(3.7) 



3.2 Dressed charges 

Given the general theory, it is natural to inquire about the form of the dressed quanti- 
ties within PN. Accordingly we would like to compute the dressed charges (the dressed 
propagators are simple and are also discussed below). 

Following our discussion at the end of subsection |2.2| on the generalization of the defi- 
nitions for dressed sub-diagrams to a multi-field field theory we define three dressed charges 
for the interaction of gravity with a compact point-like object. The (j) charge is usually 



9 Quantum Mechanically the wave-number k is equivalent to momentum, and this is how k is usually 
referred to, but in Classical field theory it is not a momentum. 
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I m 



{WnG)- 1 ^ (k • qS lm - 2k / q m / dt 



\k,i 

^> = -(SttG)- 1 ^^ • q%) / di 



\k,i 
n lm 

^)*AAA, = -(327rG)- 1 (q ■ k*ij <W - 2q ■ k<5 a fy m - (q,^ - q^)^ 
, ' q ' ^ +2q,-k,„ <5 j7 - 2q m k ; 5y + 2q m k J 5, ; - 2q,k^ jm , - 2q J k ; d~,; m ) / 



Figure 17: Static bulk vertices obtained from the expansion of the Hilbcrt-Einstein action (|3.3| ) 
and gauge fixing term ( pjj ) but restricting to time-independent terms. The undetermined wave- 
numbers flow into the vertex. 

referred to as energy, the A charge is the energy current (or alternatively, momentum dis- 
tribution) and the a charge is the stress. Together they describe the full energy-momentum 
tensor (in space-time) . In analogy with figure |3| we define the 3 dressed PN charges in figure 



19 . All the dressed charges describe the changing apparent charge due to non-linear bulk 
interaction at large, but not infinite, distances. 

In PN we have a quadratic retardation vertex for all 3 fields just like in our scalar field 
example. Since the 3 fields have different tensor characters they cannot mix, and we have 
exactly 3 dressed propagators each one defined as in figure]^ (bottom). Physically they all 
correspond to full relativistic propagators, even though A and a are spatial rather than 
space-time tensors. 

The Feynman diagrams required for the computation of the dressed energy distribution, 
p, 10 and the stress, s 13 , up to 2PN as well as the dressed momentum distribution, j l , up to 
1.5PN are shown in figure 20, figure P and fi gure El 11 The PN order was chosen as the 



one required for the 2PN effective action, apart for the stress where we chose to compute 
an additional PN order. All the results below apart from the last term of the dressed stress 
2PN were tested and confirmed against the known expression for 2PN effective action. In 
addition we shall show in the next subsections how to use the dressed charges to calculate 
diagrams beyond 2PN. 



10 In this section we shall shorten the notation and write p rather than pd r - 

11 Alternatively, one may count PN orders relative to the leading order. With this convention the dressed 
energy is computed to order +2PN relative to leading, while the dressed momentum and stress are computing 
to +1PN beyond leading. 
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s 



(167rG)- ld ^8 lm J dtdf 



d t , Im 





q, Im 



liWnG)- 1 ^|(2q^, m - S lmqi ) J dtd t 




dt, Im 



i{8ttG) ^(q/zm - feqj J dtd t 




Figure 18: Time-dependent bulk vertices obtained from the expansion of the Hubert-Einstein 



action (3.3) and gauge fixing term (3.6). Time derivative above the propagator indicates the 
direction on which it acts. The undetermined wave-numbers flow into the vertex. 



Given the Feynman rules of the previous subsection, one can write down the expressions 
for the Feynman integrals. The integrals which are essential for the evaluation of the loop 
integrals and the Fourier transforms are listed in Appendix |A|. Note also that since we 
compute the dressed sources up to a definite order in the PN expansion there is no need 
to keep the vertices of figure |l^ as they are. Instead, one has to expand each such vertex 
keeping only those powers of 'V which contribute to the PN order under consideration. 

We proceed to present the results of all the diagrams. We start from the dressed energy 
distribution. Let us denote 



A := Gm 



,r(¥)r(¥) 



d-4 



r(|) 



(3.8) 
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/ k />(t,k)0(t,-k) := « 




/di/ k j'(f,k)A:(^k) := 



J/*/ k a«(t,k)o- <i (t,-k) :=|j| 




Figure 19: The diagrammatic definition of the dressed charge distributions in PN j l and as 
the one point function for 0, A4 and er^ respectively in the presence of a single source. Retardation 
vertices are forbidden from the external leg. 




Figure 20: Feynman diagrams which contribute to the dressed energy distribution up to 2PN. 
Diagram (a) is leading (order OPN) while the rest are 2PN. 



fj%k)Mt,-k) = 




Figure 21: Feynman diagrams which contribute to the dressed momentum distribution up to 
1.5PN. Diagram (a) is leading (order 0.5PN) while the rest are +1PN (1.5PN). 



then 



-m I dt— , „ 4>(t, x(t)) 



(d-3)VT^ 



Sp(t,x) 



m I 1 + 



d- 1 



3d - 5 



2(d - 3) 8{d - 3) 



v A S(x-x(t)), 
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y kSij (t,k)ov(t,-k) 






/ >wv 



^>wv 



(a) (b) 



(c) 



(d) 



(e) 



(/) 



Figure 22: Feynman diagrams which contribute to the dressed stress tensor distribution up to 
2PN. Diagrams (a,b) are leading (order 1PN) while the rest are +1PN (2PN). 



1180(6) = A^/^ye 



Ggfflc) = A 



d-2 



dt i < 

k 



ik-x(t) \k\d-3^_^ 



(3.10) 
(3.11) 



fig.20(d) 



(lGnGym 6 (d - 3)< 



dt 



5p(t,k) 



24 (d-2)(d-5)j Jk 
(16^G) 2 m 3 (d - 3) 2 



24 



(d-2)(d-5) 



ik <^k 2 1{\k\)<j){t,-k) 
e- i% - 3 ^ P I{\k\) , 



with 



r(4-d) T(d/2-3/2) 2 r(d-3) (k 2 
(47r) d - 1 r(d/2 - l)r(3d/2 - 9/2) I T 



d-4 



(3.12) 



(3.13) 



Apparently, the above expression possesses a pole when d = 4, and thus one needs to 
introduce a counter-term which inevitably leads to a logarithmic behavior of the dressed 
energy distribution with scale. Yet this pole unphysical. Indeed, applying a Fourier trans- 
form to equation ( 3.12| ) yields 



Sp(t, r) = (8vrG) s 



m 



(d-3) 



i7-3d 



where 



d-2 



fi d _ 2 y (d-2)(d-5) 

r = x — x(t) . 



d-l 

2ir— 



(3.14) 



(3.15) 



r(¥) 

and this expression is regular for d = 4. 

Alternatively, one could stay within the wave-number space (fc-space) by introducing 
the following counter-term 

G 2 m 3 I' .f ^, , G 2 m 3 



c.t. = / dt / e-' lk ^ k 2 cj)(t, -k) 

6(d - 4) 7 7 k n ; 6(d - 4) 



dtn<j)(t,x(t)) . (3.16) 



However, as explained in j|] such a term can be removed from the Lagrangian by an 
appropriate field redefinition as it is proportional to the leading order equation of motion 
for the NRG field d> and hence it is a redundant term. 12 



12 In quantum field theory such an objects is called a redundant operator referring to its action on the 
Hilbert space of states. In the classical theory it is not an operator and calling it a term is more appropriate. 
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The contribution to the energy distribution of the diagrams which contain time deriva- 
tives is given by 



figMe) 



- I dt I < 



-ik-x(t) 



|d-5 



k 2 v 2 - {3d -5)(k- vf - 4(d - 2)(ik ■ a) <f>(t, -k) 



d-2 



A dt e 



-ik-x(t) I jud— 5 



ik-a + 2{k-v) 2 (j)(t,-k) (3.17) 



Combining altogether we get up to 2PN 



p(t, k) e ik ^ = m I 1 + 



d- 1 



3d -5 



2(d-3) 8{d-3) 
Ad 2 - 17 d + 19 



(16vrG) 2 m 3 (d - 3) 2 



24 



(d-2)(d-5) 



fc 2 /(|fc|) 



d- 3,,- 



x((/ _. :(:i X\k\^v 2 + ^\\kr\k.vY. (3.18) 

where the exponential on the left hand side really belongs on the right hand side and was 
moved from there to achieve a "cleaner" form. Transforming back to coordinate space 
yields 

/ ( d—l-y 3d — 5 4 \ 
p(t, r) = m I 1 + — — — v + — r- — —v o(r) 



2{d-3) 8{d-3) 

2 



fid- 



+ (8vrG) J 



d-2 

! 'i 

d-3 



4-2d 



(d-3) 



f\ 7 - 3d . (3.19) 



,O d _ 2 y (d-2)(d-5) 
The results for the dressed momentum distribution up to 1.5PN are 

m v 1 



fig©' 



df 



Ai(t,x(t)) 



1 



5f(t,x) =m« 1 1 + y S(x-x(t)), 



Sf(t,k) = mv l ( 1 + -v 2 ) exp(-zfc • , 



fig.|3(6) = -A 
figH(c) = A 



d-2 
2 

d-3 



df 



; -ifc-2(*) |fc| d - 3 {/ -A(t,-k) 



8(d- 2) 



dt / e - i ^W|A?| d - 5 (fc 2 ^ + (d-3)(A?-t;)A; i )v4 i (t,-fc) 



(3.20) 



Altogether we obtain up to 1.5PN 



f(t,k)e ik - s ® =mv i [l + -v 



1 



4d 2 - 17 d + 19 
8(d- 2) 



\d-3„ 



+ ^^Al^Cfc-tO**. (3.21) 
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In coordinate space the dressed momentum distribution is 

2 



(v ■ r)r — 2v 



\r\ i ~ 2d . (3.22) 



d- 2 

Finally, the results for the dressed stress charge up to 2PN are given by 



fig-! 



*.> = ?/ 



dt 



V V 



VT^ 2 



(Tij(t,x) 



ds ij \t,x) = mv i v j 1 + -v 2 5(x -x(t)) , 



5s lJ (t, k) = m v % v 3 ( 1 + -v 2 ) exp(— ik ■ x(t)) , 



(3.23) 



fig| 



x (d-3) 2 f , / d-1 ,, 
|(6) = X-^—J-^ I dt ( 1 + - — r u 2 



16(d - 2) 



d-3 



|A:| d - 5 kikj-k^Sij a ij (t,-k) 



(3.24) 



fig#) = A 



32(d-2)J y k 1 1 

(3d-5) 2 . „ 2 , r2 ^ 2 , 8(d-l)(d-2) 



d-3 



-(fc-^ + A; 2 ^ + 



d-3 



ifc • d 



5ija ij {t,-k) (3.25) 



figj2g(d) = A 



d-2 



dt / e - **"^ |jfe| d - 5 
2 ■/ Jk 



~{k-v) 5ij + (ik ■ a)5ij - ikittj 



a ij (t,-k) 

(3.26) 



fig|2|(e) = - j dt J e -i*<t)\k\d-5 J( d _ 3)^2^2^. _ _ 2 ( d _ 2)fc 2 ^-] ^'(t, -jfe) 

(3.27) 



fig-i(/) 



.^^| y dt y e -**(t) |jfe|d-5 

^P* 2 3)^-4) ^ . 2 _ d 2 -7d + n fk 

8 8 v 1 4 v 7 

.^^| y dt y ,,-*•*(*) |jfe|d-5 



8ijU io (t, —k) 



d a ^ v 2 kikj + — — — —(A; • v) 2 kikj - \k 2 ViVj 
o o 4 



t'J 



Mzl j dt J e-**^*) |fc| d - 5 



^ + (d-3)(d-5) (jg 



<7« 

(3.28) 
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As a result, the dressed stress charge up to 2PN is given by 

k) =mv i v> (l + \v 2 ^ - ^^[^ 19 + yBV'j |fe| d " 5 A 

8(rf - 2) V 2 2 v ; v fc2 i 11 



+ 



'X8 13 



(d-1 )(d 2 + 8d-21 ), r . a d 3 + 2d 2 - Yld + 7 r 

- —. r -(k ■ vY H — r ia ■ k 

16(d-2) y ' 8(d-2) 

(d_ 3 )2 ^_^\^ |d _3 A ^ ; (329) 



8(d-2) V 2 

where (ij) denotes symmetrization with respect to indices i and j with factor 1/2 included. 
In the coordinate space we obtain 

.«(t,f) = mtV (l + i» 2 ) 8(f) - (<j i 3 ^/l 3) /(Of"a" - 2/(r)vV 

- 4f/-") /(r)( " 2 - 2(rf - 2)( 7 )2)i ' J 

- ^ — -/(r) ( (d - l)(u ■ f) 2 f *f J ' - 2(v • f)v*f j + - ( 1 - — J 5 ij - r i r j + (a ■ f)r i r j 
d 2 V 2 V 2 / 

+ 2(l + 2K7-S) ^P-0^- (3-30) 

with 

(«£)'■ <-) 

3.3 Skeletons for 2PN and 3PN 

We shall now demonstrate how all the bare diagrams of the two-body effective action up 
to 2PN, order by order transform into their dressed form, including their corresponding 
skeletons. 

At OPN (Newtonian order) a single diagram (figure |l|) contributes representing the 
interaction of two tree-level masses. 



At 1PN (Einstein-Infeld-Hoffmann Lagrangian) the four diagrams shown in figure 23 
contribute. The first is a v 2 component of the mass vertex and can be considered to be a 
trivial dressing of the energy. The second is the tree-level interaction of two currents. The 
third represents a propagator dressing (retardation) of the Newtonian potential. Finally 
the fourth is due to a non-linear world-line vertex which accounts for the gravitating nature 
of potential energy. Altogether none is irreducible with a topology other than Newtonian. 
We note that the last diagram is of order O (G 2 m 3 u°) where m represents a typical mass 
and v a typical velocity while all the previous diagrams are of order O (G 



m v 



At 2PN Gilmore and Ross []9| found 21 diagrams. We shall see that many can be in- 
terpreted to represent dressing effects while only one is both non-factorizable and dressing- 
irreducible. Indeed || find its computation to be the core or essential computation at this 
order. 
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Figure 23: The four diagrams which contribute to the two-body effective action at 1PN 
Einstcin-Infcld- Hoffmann Lagrangian. 



The 



-X- 



-x- 



(a) 



(b) 



(0 



(d) 



wvw\ 



(/) 



Figure 24: 2PN Diagrams contributing at order O (Gm v 1 ) (following M) 



The six diagrams (a-f) shown in figure 24 contribute at order O (Gm?v 1 ). Diagrams 
(b,c,e) represent propagator dressing to diagrams from lower orders, while the rest involve 
only tree-level v-dependent vertices. 

At order O (G 2 m 3 ?; 2 ) there are 10 diagrams shown in figure 25. Diagrams (a,b,c) 
are V-shaped and as such factorizable and at least one factor is the Newtonian potential. 
Diagrams (d-j) are Y-shaped and as such are dressing-reducible. Four diagrams (d,f,g,i) 
represent mass dressing, the two (e,h) represent current dressing, while finally (j) can be 
thought to represent stress dressing. 

At order O (G 3 m 4 u ) there are 5 diagrams. Both (a) and (c) factorize into 3 Newtonian- 
potential factors. Diagrams (b) represents a mass dressing (the circled piece) while (e) 
includes two a dressing sub-diagrams. Finally diagram (d) is the one and only truly irre- 
ducible diagram at 2PN. Actually for some yet-unexplained reason the computation reduces 
after several steps to a square of the master one-loop integral. We speculate that this is 
special to the GR action (and is not generic to classical field theories) and in particular to 
the gauge symmetry which may relate this diagram to other diagrams at 2PN. 

As a result of this analysis we can extract all the non-factorizable skeletons up to 2PN. 



These are listed in figure 27, where the skeletons are labeled according to the PN order 
in which they first appear. By evaluating the dressed diagrams we successfully tested 
our expressions for the dressed charges from the previous subsection against the known 
expressions for the effective action. 

3PN. As a step towards the determination of 3PN we extend the list of skeletons up 
to that order using the classification of possible topologies in figure [15| and some knowledge 
on the Feynman rules of PN. Taken together with the evaluation of the dressed vertices 
(partially obtained in the last subsection) this list of skeletons leads the way to a determi- 
nation of the 3PN part of the 2-body effective action. We note that the 3-loop topologies of 
figure [l5| are not realized at 3PN. The reason is the absence in PN of certain bulk vertices 
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Figure 25: 2PN Diagrams contributing at order O (G 2 m 3 u 2 ) (following §). 






Figure 26: 2PN Diagrams contributing at order O (G 3 m 4 v° ) (following §) 



such as </> 3 . 

3.4 Computing beyond 2PN 

In this subsection we explicitly demonstrate the economizing ability of proposed approach 
in several cases 

• A certain economization during 2PN calculation. 

• Using the dressed energy distribution computed above in subsection |3^ to calculate 
terms of the Newtonian interaction type of order 3PN and 4PN. 

• Similarly, using the computed dressed momentum distributions we calculate a current- 
current interaction term of order 3PN. 

• Finally, we use our highest order results for the dressed stress in order to compute 
3PN terms by attaching the external a leg to the second compact object. 

An application at 2PN. The first non-trivial example appears during the computa- 
tion of the 2-body effective Lagrangian at order 2PN. Indeed the dressed stress sub-diagram 
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Figure 27: A listing of all non-factorizable skeletons appearing in the dressed Post- Newtonian 
perturbation theory up to 3PN, listed by the PN order at which they first appear. 




Figure 28: The diagrammatic representation of the recursive integral equation satisfied by the 
dressed energy. 



in figure ^(b) appears twice at 2PN - both at figure |25|(j) and at figure ^6j(b) , and of course 
it is enough to evaluate it once. Another point of view is to obtain both dressed energies 
corresponding to figures |20|(c,d) from the recursive relation shown in figure [2^ after sub- 

|a,b). 



stituting in the leading contribution to the dressed stress from figures 

Proceeding beyond 2PN we shall now compute certain 3PN and 4PN terms in the 2- 
body effective action building on the results of the previous section. Pictorially these terms 
are shown on figure 29, where the bubbles represent the dressed charges, and we explicitly 
indicate on each bubble which piece of the dressed charge is essential to the computation. 
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(a) 



(b) 



(c) 



4PN: 





(d) 



Figure 29: The diagrammatic representation of terms beyond 2PN in the 2-body effective action: 
(a) a 3PN p-p term, (b) a 3PN j-j term, (c) a 3PN s-s term and finally (d) a 4PN p-p term. 



Building on the Feynman rules of the previous section, we obtain 



fig.H(o, d) = 8ttG 



d-3 
d-2 



dt 



fig.|29|(6) = -16vrG / dt 



p 2 (t,-k)pi(t,k) 
k k 2 
j 2 (t,-k) • ji(t,k) 



k k 2 



fig. 29 



8irG 



S dt L 



s l 2 J (t,-k)P ijM s? l (t,k) 



(3.32) 
(3.33) 
(3.34) 



The 3PN j-j and s-s terms. Substituting (|3T2l|) , fl3~29p and d = 4 into (|3.33| ),(|3~3l 

yields 



fig-29(6) 



Gm 1 m 2 _ 2 2 

dt [Vi ■ V2) V-y V 2 

r 

Gm x m 2 ( Gm 2 2 Gm x 2 N , , - w . -.n 

ai — ; I v x H u 2 [ l\v\ ■V2) + (v\- r) (v 2 ■ r)\ 



2r 



+ 



7 li 
32 



dt 



(vi -v 2 - 3(«i -f)(v 2 ■?)) , (3.35) 
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fig.H(c) = / dt^f^ [(* • v 2 f - v 2 v 2 ] v\ + Jdt 



/ -> -* \ 2 203 9 9 
2^1 "°2) + "]J U 2 U 1 



+ / dt 



dt 



G mf m 2 



~(^2 • r){ai ■ v 2 ) - -^-vlifli ■ r) - -(ai • r)(v 2 ■ r) 2 



G mf m 2 



+ 



245 2 _ „ 2 1 
g-^(« 1 .r) 2 + i 

1 f G 2 m\ m 2 



^-^v^ivi ■ r) 2 + ^vf(v 2 • r) 2 + -(v\- r) 2 (v 2 ■ r) 2 + -^(vi ■ v 2 ^ 2 

o 



1 

12' 



eft- 



[(«! • U2)(ui • ?){v 2 ■ ?)] + (1 <-> 2) . 



(3.36) 



where f(t) = x 2 (t) — x\(t) denotes the radius vector between the particles, and we have ne- 
glected the contact term proportional to 5(r) since in CLEFT particles are widely separated 
from each other. 

The p-p terms. Let us now evaluate figure [291(b). First note, that since equation 
( 3.12 ) contains a simple pole in d = 4 one concludes that the integrand of ( |3.32|) possesses 
a double pole, and thus in contrast to the previous case ( 3.14| ) a Fourier transform will not 
lead to a regular expression. 

Therefore, we expand equation ( |3.12j ) around d = 4 



fig.[20](d) = G 2 m 3 r L~ e j dt J e~ ik ' s ® k 2 <f>(t, —k) 



~ + A (- 1 + 2 ^ " 21n(87r 2 ) + 21n(£ 2 L 2 ) ) + 0(e) 
be 12 \ 



(3.37) 



where 7 is Euler's constant and motivated by the QFT renormalization approach we intro- 
duced an arbitrary length scale L and the following definitions 13 



e = 4 — d , Gm = Gm r L e . 



(3.38) 



In order to eliminate the divergence as e — > 0, we add the following counter-term to 
the effective action 



c.t. ■ r j ,11 A< >(/.. r(/)) -( j dt j e~ ik - S{t) k 2 c/)(t, -k) 



(3.39) 



where 



c = L e ( c r 



G 2 m 3 r 
6e 



(3.40) 



13 Index "r" stands for " renormalized" , though in the current work we will not encounter divergences 
associated with non-trivial RG flow of the mass. 
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Since the renormalization scale L introduced above is arbitrary, we must have 

dm r " 



= L^ = L- 
dL 

dm r 
L — = e m r 



-e m r + L- 



dL 



= L 
=> L 



dL 
dc 
IE ' 
dc r 
~dl~ 



-eL~ 



ec r + 



G 2 mf\ / dc r G 2 m 2 dm r 

G 2 ml 



(3.41) 



Apparently, the theory exhibits a non-trivial RG flow. But as argued in |B[ this scaling 
is not physical in nature and can be removed by a suitable field redefinition which is 
tantamount to a coordinate transformation. Indeed, combining ( p.37 ) with ( |3.39| ) yields 

(;2l "' : ' -l + 2 7 -21n(87r 2 ) + 21n(£ 2 L^ 



Sp(t,k) 



Cr 



12 



-ik-x{t) p _ ^ 3 42 ) 



In coordinate space this expression is given by 

G 2 m 3 



6p(t,f) 



12 



(-l + 2 7 -21n(87r 2 )) 



6(7 



1 G 2 ml 
2tt r 5 



where f = x — x(t) and the relation 



d 



lnfe 2 = lim (k' 
a^o da 



(3.43) 



(3.44) 



was used in order to apply the Fourier transform formulas of Appendix |A|. 

However, the contact term in the above expression for 8p can be eliminated by the 
following field redefinition 

' -l + 2 7 -21n(8vr 2 )) 



+ 4ttG 



12 



6(r). 



(3.45) 



As a result, we reproduced equation ( p. 14 ) for d = 4 and also obtained the regularized 
expression for p(t, k) in four dimensions 14 



p{t,k)e ik - s{t) =m(l + lv 2 + lv 4 )+lGm 2 \k\ (l5v 



7T 



G 2 m 3 ? 



2 T 2\ 



k z ln(k z L 

(3.46) 

A comment should be made regarding the field redefinition ( 3.45| ). Obviously, such 
a shift introduces an extra term to the world- line action as a side effect. This term is 
proportional to the second derivative of the NRG field (f) with respect to time. However, it 
can be removed by including additional counter-term given by 



c.t. 



dt -^(t,x{t)). 



(3.47) 



We omit index "r" to avoid abuse of notation. 
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We will not elaborate the details of this elimination as they are irrelevant to the computa- 
tion of the 4PN term. 

Substituting ( 3.46 ) into ( 3.32|) leads to 



fig||(a) = ^ 



and 



Gmim 2 2 4 
at V\V2 



+ -ldt 



G 2 m\ m\ 



2^2 + 



(v 2 -r) 2 lGm 2 



3 r 



+ (1 



(3.48) 



fig®0 = | 



1 

+ 3 

1 

+ 3 



+ 



7T 



(ft 



Gmim.2 



44 + ~ 



+ 



G 2 m| mi 4 
* ^ v i 

mf m2 4 
5 v 2 



^ 2 + ^ 
2 2 2 



2 ^ lGm 2 



7 2+ (^)_. + 
2 1 2 3 



3 r 
1 Gmi 



, G 3 m? m 2 , r/ _ _ N o 



59 2 2 , 87 



■•- 1 -r) 2 



15 



+ it (^2 • r) (vi ■ r) - 6(v 2 ■ r)(vi ■ r)(vi ■ v 2 )} 



dt \ — 2 - 



2 /23 
vf — 



dt 



G 4 m\ m\ 



Hi 

23 



87 - 8m(r/L) + (v x -r 



87 

2 ' 



hx(r/L) + (v 2 ■ ?f 



- - + 2 7 + 21n(r/L) 
-- + 2 7 + 21n(r/L) 



if 



dt 



G 5 



1 "" 2 [3-2 7 -21n(r/L)] . 

(3.49) 



To derive this expression we first applied ( 3.44| ) and then used transform Fourier master 
integrals of Appendix [A|. Appearance of an arbitrary renormalization scale L in the above 
expression should not be taken as a strong evidence for a non-trivial RG flow, since if other 
terms in the full 4PN potential are considered then cancelation of logs might occur. 
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A. Master integrals 



In this appendix we present master integrals which we found useful during the calculations 
presented in the text. We start with 



J 



d d q 



( k 2)d/2-a-/3 Y{a + (3- d/2) T(d/2 - a)T(d/2 - 0) 



(27r) d (q 2 ) Q [(q-k) 2 ]' 3 (4vr)<V 2 T(a)Y(p) 

d d q qi d/2 — a 



T(d-a-(3) 



J; 



( 2vr )d( q 2)a[( q _ k )2]/3 d-a-f3 



(A.l) 
(A.2) 



J; 



d d q 



(2n) d (q 2 )«[(q-k)Y 

1 r(a + p - d/2 - 1) T(d/2 - a + l)r(d/2 - /?) 



(47r) rf / 2 r(a)r(/3) 



T(d-a- P + 2) 



{d/2 -a + l)(a + p-d/2- l)hkj + (d/2 - P)—6ij 



(k 2jd/2 



-a-P 



(A.3) 



J, 



d d q 



% ?j g& 



M// ; (27r) d (q 2 ) Q [(q-k) 2 ]/ 3 

T ( a + p_ d / 2 _ i) r(d/2 - a + 2)r(d/2 - /?) 



( 4vr )rf/2 



r(«)r(/?) 



r(d-a-/9 + 3) 



{d/2 - a + 2) (a + 0-d/2- l)kikjk k + (d/2 - P) — (5ijk k + + ^Ay) 



(A.4) 



ijkl 



d d q 



Qi Qj Qk Ql 



r( d / 2 - a + 2)r(d/2 - /3) 



(27r) d (q 2 ) a [(q-k) 2 ]^ 
r(a + P - d/2 - 2) 



( 4vr )rf/2 



r(d - a - P + 4) 



.2\2 



r(a)r(/3) 



[(d/2 - P)(d/2 -p+ l)(Sij8 k i + 5 ik 5 jt + 6uS jk )- 



+ (Sijk k ki + 5 ik kjki + 5ukjk k + 5j k kiki + 5jihk k + Skihkj) — 

x (a + P - d/2 - 2) (d/2 - a + 2)(d/2 - /3) 
+ kikjk k k x (a + P - d/2 - 2)(a + P - d/2 - l)(d/2 - a + 2) (d/2 - a + 3)] (A.5) 

As noticed in || the above integrals ( |A.2| )-( [A5| ) of vector nature can be reduced to a scalar 
integral ( |A1| ) on the basis of their transformation properties under rotations. Yet such a 
reduction becomes involved when the rank of the tensor under consideration increases, and 
therefore we choose to list all those which were relevant to this work. 

In order to evaluate the above integrals we proceed as follows. We first apply the 



generalized Feynman parametrization [4C] 

N i N 

ni=in^(E 



i=l 



Xi — 1 



nu—l 



r(mi + m 2 + ... + tun) 



[y>iA] Em< r(mi) • • • r(mjv) 



(A.6) 



- 33 - 



with N = 2 and m\ = a, m<i = (3. Next we integrate over one of the Feynman parameters, 
e.g. X2, and subsequently redefine the undetermined wave-number q — > q+ x±k. 
Now, in order to integrate over q, we build on the following formula 

f A 1 _ z ~ d ' 2 nn - d/2) d/2 _ n 

J (2vr) d (zq 2 + A) n (4vr) d / 2 T(n) ' { ' 

computed by means of dimensional regularization. Thus, for instance, differentiating it 
with respect to z a definite number of times and setting z = 1, one obtains 

f A q 2 _ d/2 T(n - d/2 - 1) d/2 _ n+1 
J (2n) d (q 2 + A)™ (4vr) d / 2 T{n) 

f d\ (q 2 ) 2 d(d + 2) T(n-d/2-2) d/2 _ n+2 

J (2ir) d (q 2 + A)™ 4(4vr)<*/ 2 T(n) ' [ } 

As a final step we integrate over x\. 

Another set of useful identities is related to the following d— dimensional Fourier trans- 
form 



d d k e ikr 1 V{d/2-a) fr 2 \ a d/2 



(2vr) d (k 2 ) Q (4vr) d / 2 T(q) V4 
Differentiating it with respect to r yields 

d'^ k, gikr= . x T(d/2-a+l) fr^^ ' 



(A.9) 
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(2vr) d (k 2 )° " * 2(4vr) d / 2 r(a) \A 

d d k kikj ikr _ F(d/2 -a + l) fSij_ + _ _ 1 x^j \ f r 2 ^ °- d / 2 - ] 



(2vr) d (k 2 )" (4vr) d / 2 r(a) V 2 v ' r 2 

tf*k kikjk, ikr _ iT(d/2 - a + 2) / r 2 \ a ~ d/2 ~ 3 



(A.10) 
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(2n) d (k 2 )° 16(4vr) a! / 2 r(a) V 4 

x [r 2 (6uxj + SjiXi + Sijxi) - (d - 2a + A)x i x j xi] , (A.ll) 
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